Non-magnetic B-site Impurities Induce Ferromagnetic Tendencies in CE Manganites 
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Using a two-orbital model and Monte Carlo simulations, we investigate the effect of nonmag- 
netic B-site substitution on half-doped CE-type manganites. The lattice defects induced by this 
substitution destabilize the CE phase, which transforms into (1) the ferromagnetic (FM) metallic 
competing state, or (2) a regime with short-range FM clusters, or (3) a spin-glass state, depend- 
ing on couplings and on the valence of the B-site substitution. While a C-type antiferromagnetic 
state is usually associated with an average eg charge density less than 0.5, the nonmagnetic B-site 
substitution that lowers the eg charge density is still found to enhance the FM tendency in our 
simulations. The present calculations are in qualitative agreement with experiments and provide 
a rationalization for the complex role of nonmagnetic B-site substitution in modulating the phase 
transitions in manganites. 

PACS numbers: 75.40.Mg, 75.47.Lx, 75.47.Gk 
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I. INTRODUCTION 

Rare-earth manganites of the form i?i_a;Aa;Mn03 
(where R (A) is a rare-earth (alkahne-earth) element) 
are typical representatives of complex oxides with 
multi-orbital interactions and a strong competition be- 
tween spin, charge, orbital, and phononic degrees of 
freedom In these materials the main competition is 
between the delocalization tendency of the eg electrons 
and the localization effects caused by the antiferromag- 
netic (AFM) coupling between the Mn t2g spins as well 
as the Jahn- Teller effects jii^i^i^ The delicate balance be- 
tween these competing tendencies produces a very rich 
phase diagram. Intrinsic or external perturbations that 
naively may seem "weak" , such as small variations in the 
carrier density, pressure, magnetic fields, and quenched 
disorder, can nevertheless induce nonlinear effects, in- 
cluding phase transitions due to the close proximity in 
energy of the competing states. This high sensitivity to 
perturbations is clear in manganites with composition 
X = 0.5, where several competing ground states with 
quite different properties, such as ferromagnetic as well 
as CE and A-type AFM states, have been identified i^i^i^ 
Among them, the CE state is known to appear in sev- 
eral narrow-bandwidth manganites. This complex state 
is formed by zigzag FM chains with AFM inter-chain cou- 
pling, and its stabilization is usually accompanied by a 
checkerboard pattern of charge and orbital order, which 
further stabilizes the CE spin structure.—i&i^ 

Metal-insulator transitions (MIT) are familiar phe- 
nomena in condensed matter physics. For a nor- 
mal metal, the addition of quenched disorder leads 
to the trapping of mobile carriers and transforms 
a metal into an insulator. However, for mangan- 
ites the reversed process usually occurs, in situations 
where quenched disorder originates from chemical sub- 
stitution or intrinsic defectsi^ Usually, two sources of 



qucnchcd-disorder effects in manganites are considered. 
One is primarily caused by the A-site disorder (al- 
loy randomness) fiSiiii while the other is induced by 
B(Mn)-site substitution^asiiOSieaiisag^^^ Theo- 
retically, the A-site disorder effects have been exten- 
sively studied in the past few years. Bond disor- 
der and on-site potential disorder were jointly or sep- 
arately introduced to model various A-site disorder 
sources^^>26^>28^^^^^^^ The quantitative 
calculations indicate that the FM metal phase is not 
much affected by the A-site disorder, while the CE/CO 
insulator easily collapses into a glassy-like state. 

Conceptually, and differently from the A-site disor- 
der, the B-site substitution induces disorder locally di- 
rectly into the Mn-0 network and it may modulate 
quantities such as n, the Cg electron-density. There- 
fore, the B-site substitution can have stronger impact 
on the physical properties of the material than the A- 
site disorder. Experiments show that a few percent B- 
site substitution, such as chemical substitution of Mn 
by Cr/Al/Ga/Ru in the half-doped CE state of man- 
ganites, will favor a phase separated (PS) state with 
FM tendencies, although these ions are different in their 
electronic structure and have different magnetic coupling 
with the Mn ionsJ^iiiii^J^iiiii^ii^i^iSi For the nonmag- 
netic substitution, for instance, 2.5% Al substitution 
in Pro.sCao.sMnOs is sufficient to convert the charge- 
ordered (CO) CE insulator into a state with FM metallic 
characteristicSi-*^ These experimental results are very sur- 
prising since (1) the substituting ions Al/Ga are nonmag- 
netic, and (2) the Al/Ga doping leads to the reduction 
of the Cg electron density, both of which are disadvanta- 
geous for the FM tendency^^ In fact, such a reduction 
of the Cg electron density is expected to favor an AFM 
insulating statci^ 

Very recently, the disorder effects by B-site substitu- 
tion were also investigated within the context of the two- 
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orbital double-exchange (DE) mode L^^i^^ The observed 
collapse of the CE/CO phase into a FM phase was ex- 
plained in terms of a density-driven phase-separation. 
The main idea is that the impurity, having a -1-4 valence, 
transfers extra electrons to the remaining Mn sites.'^^ 
Thus, in this context it is straightforward to understand 
the FM tendency since Cg densities larger than half, 
nr>0.5, usually are associated with FM phases in man- 
ganites. However, this density-driven phase-separation 
idea can not explain the experimental fact that trivalent 
substitutions, such as Cr/Al/Ga, are expected to reduce 
the eg electron density on the remaining Mn-site (ni.<0.5) 
rather than increasing this density. This reduction of the 
Cg density is quite nontrivial since the electronic density 
is one of the most important factors to determine the 
ground state, especially around the half-doping region. 
In fact, the previous theoretical investigation predicted 
that Al/Ga would not lead to a FM tendency for half- 
doped manganites based on the previously described idea 
of a density-driven phase-separationi^ Therefore, it is 
necessary to re-investigate the B-site-substitution disor- 
der effects in half-doped manganites to understand these 
puzzling experimental results. 

In the following, we will study the effects of B-site sub- 
stitution on the stability of the CE/CO state. Since even 
Al/Ga substitutions can induce the FM tendency, it is 
reasonable to consider first a non-magnetic impurity for 
simplicity. The main results found in this paper are that 
the B-site disorder can induce two main effects: (1) The 
lattice defects by B-site substitution break the CE-zigzag 
chains and frustrate the charge ordering which destabi- 
lizes the CE phase and induce the FM tendency (the FM 
state is close in energy to the CE phase in half-doped 
manganites); (2) The reduction from half-doping of the 
Cg electron density suppresses both the long range FM 
and CE tendencies, since it prefers the C-AFM state. 
These two competing effects can trigger a phase tran- 
sition from the CE/CO phase either into a short-range 
FM cluster state with relatively strong FM tendencies, or 
into a spin glass state. Due to these competing effects, 
an optimal B-site substitution level for the stabilization 
of the FM tendency is found to exist. 



II. MODEL 

In our investigations, we consider a two-orbital model 
defined on a two-dimensional Lx L square lattice (L — 8, 
and using periodic boundary conditions) with a Hamil- 
tonian given by, 

H = - ^ iapAaadjPa + J AF ^ • Sj 
<ij>,a,fl,a <ij> 

-Ja ^ • Sj + X^^iQupi + Q2iTxi + QsiTzi) 

i i 
i i 



where the first term is the two-orbital DE interaction, 
^Lo- (diaa) is the Creation (annihilation) operator for 
an Cg electron with spin a in the orbital a {(1^2 
or d^^2_^2) at site i. The hopping amplitudes be- 
tween nearest-neighbor (NN) sites < ij > are given by 
tla = -V3t:, = -V3tL = ittt = 1 for V = X, and 
tL = V3tl, = V^tl = 3tl ^IfoTv = y. The sec- 
ond term is the AFM super-exchange (SE) interaction 
between the NN t2g spins S. In the third term, the Hund 
coupling Jjj (> 0) links the Cg electrons with the 
spins S (assumed classical and normalized as |S| = 1). 
For simplicity, we consider here the Hund coupling in 
the widely-used limit of Jjj — > oo. The fourth term is the 
electron-phonon coupling, where A is the dimensionless 
coupling constant, Q are the phononic modes {Qi is for 
the breathing mode, Q2 and Q3 are for the Jahn- Teller 
modes), and r is the orbital pseudospin operator. The 
fifth term in the Hamiltonian is the elastic energy of the 
phonons.'^* For simplicity, and as in many other previous 
investigations, the phonons here will actually be consid- 
ered just as classical lattice distortions, /i in the last term 
is the chemical potential to tune the eg electron's density. 

To introduce a non-magnetic impurity (without 3d 
electrons) as the B-site substitution, we will assume that 
the impurity has no contribution to the electron conduc- 
tivity and, thus, consider it as a lattice defect. Thus, 
the DE, SE, and Jahn- Teller couplings around the im- 
purity can be ignored, retaining only the elastic energy 
of local phonons. These localized defects distinguish the 
B-site substitution models^i^ from the A-site disorder 
models24i25^S^S^^^^^^ ^hich the disorder 

effects were applied to all sites. Therefore, the topological 
structure of the lattice defects, which is absent in A-site 
disorder cases, is especially important in patterning the 
electron configuration of the original CE/CO state. 

Our model Hamiltonian Eq. (1) is studied via a com- 
bination of exact diagonalization and Monte Carlo (MC) 
techniques: classical t2g spins and phonons evolve fol- 
lowing the MC procedure; and at each MC step, the 
fermionic sector of the Hamiltonian is numerically ex- 
actly diagonalized. The first 10^ MC steps are used for 
thermal equilibrium and another 10"^ MC steps are used 
for measurements. More details about this widely-used 
two-orbital Hamiltonian and the MC algorithm can be 
found in Ref. [3- In the present calculation, first we con- 
sidered averages over several defect configurations. How- 
ever, we observed that for a dilute distribution of defects, 
namely with defects not in close proximity to each other, 
the results of the calculations are almost the same for 
different configurations. Therefore, here only two de- 
fect configurations were used for each parameter point 
in most simulations, except for the cases of the phase 
diagram and density-of-states (DOS) for which only one 
configuration was used. In addition, since defects should 
have the same probability to occupy the two types of 
sites in the CE phase, the bridge sites (Bl) and cor- 
ner sites (B2),'^^ the same number of defects on these 
Bl and B2 sites are arranged in our 8x8 lattice. The 
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FIG. 1: (Color online) (a) Jaf-A phase diagram at 1=0.5 in 
the absence of B-site substitution (clean limit). Here, point 
A refers to the particular set of couphngs Jaf=0.1 and A=1.4 
that is emphasized in our analysis, (b) The eg charge density 
vs. chemical potential /i, at Jaf=0.1 and A=1.4 (point A in 
Fig.l (a)). 

8x8 lattice is enough to describe the prominent phenom- 
ena of the B-site substitution, such as the destabiliza- 
tion of the CE order, the phase separation, and the FM 
tendency, since the lattice size effects are mild in these 
phenomenai^i^^ All the simulations are performed at a 
low-temperature fixed at T=0.01, which is low enough to 
describe ground state properties. To characterize differ- 
ent spin orders, the spin structure factors are calculated 
by performing Fourier transforms of the real-space cor- 
relation functions4S 



III. RESULTS 

A. Clean limit phase diagram 

First, we will briefly review the phase diagram of the 
model used here when the A-site doping is a; = 0.5, and 
in the clean limit. The result is shown in Fig. 1(a). This 
phase diagram in the (Jaf, A) plane can be divided into 
three main regimes. When both Jaf and A are small, 
the DE interaction dominates and favors the FM metallic 
state (regime denoted by "FM"). With increasing Jaf, 
i.e. enhanced SE interaction, an appropriate combina- 
tion of Jaf and A leads to the CE/CO insulating state 
(regime "CE/CO"). Moreover, there exists a parameter- 
space region with coexisting FM order and charge order 
(regime "FM/CO") where Jaf is small and A is large. 
This phase diagram was established before and verified 
experimentally, thus the reader is referred to previous lit- 
erature for more details.-'-'- Here our attention will con- 
centrate on the CE/CO regime near the boundary with 
the FM metallic phase. Considering point A (Jaf = 0.1 
and A = 1.4) as an example (in the rest of the paper, Jaf 
is fixed to 0.1 unless otherwise stated), Fig. 1(b) presents 
the ground state density in the vicinity of n = 0.5, as a 
function of the chemical potential fi. A clear plateau with 
n — 0.5 indicates a fairly stable CE phase. The other two 
plateaus correspond to two other phases: the FM state 
at n > 0.5 and the C-type AFM state at n < 0.5, indi- 



cating the importance of the charge-density variation in 
driving the phase transition. The transitions between the 
three phases are abrupt when varying the chemical po- 
tential /i, suggesting density-variation-driven first-order 
phase transitions, at least in the small clusters we have 
studied in this effort ,^i^i^i22i36^ 



B. Effect of lattice defects 

Now let us investigate the effect of B-site non- 
magnetic substitution in manganites of the form 
-Ro.5^o.5Mni_j,i3j^03. Such a substitution will lead to 
the appearance of lattice defects and, simultaneously, a 
variation of the Cg electron density. To clarify their re- 
spective roles, here we first address the effect of the lat- 
tice defects. For such purpose, the substituting cations 
are assumed to be -1-3.5 in charge to keep the average Cg 
charge density for the remaining Mn-sites (n^) invariant, 
i.e. rir = 0.5. 

We have observed that for a given appropriate substi- 
tution level y, the CE/CO state will turn into a state 
with strong FM tendency. Fig. 2(a) shows several typ- 
ical spin structure factors S'(q) evolving with y at the 
point (Jaf=0.1, A=1.4). The FM order at q=(0,0) 
emerges a.t y > 0.031 and it is enhanced quickly up to 
S'(q)--0.12. On the other hand, the E-type AFM order 
at q=:(7r/2, 7r/2) and C-type AFM order at q=(0, tt) are 
rapidly suppressed when y reaches 0.063, indicating that 
the CE spin order is destroyed by the lattice defects. This 
result is qualitatively similar to the result in Fig. 3(d) of 
Ref. [3^ However, the origin of this CE to FM transition 
is by lattice defects (to be explained below), instead of 
the Cg density enhancement proposed in Ref. [3a. The CE 
destruction can be further understood by observing the 
MC snapshot of a spin configuration at y—0.09A substi- 
tution. As shown in Fig. 2(b), there is no trace of any CE 
chains, the CE phase is converted into a state consisting 
of small FM clusters with various orientations, similar to 
the results in previous A-site disorder effortsi^ 

The fundamental reason for the CE/CO destabiliza- 
tion can be understood based on the breaking of the 
zigzag FM chains and concomitant charge frustration. 
Considering first the magnetic order, the CE phase con- 
sists of zigzag FM chains that are easily cut down by 
lattice defects, leading to a substantial increase in the 
kinetic energy. However, the competing FM phase has 
a two-dimensional (three-dimensional in real case) char- 
acter, which is much more robust against lattice de- 
fects. For a pure system the Jahn- Teller coupling fa- 
vors the long-range staggered CO pattern. In real cases, 
the B-site substitution should be randomly distributed 
between the Bl sites (with higher charge density) and 
B2 sites (with lower charge density with the same prob- 
abilities. This randomness of B-site substitution can 
break the original CO state, causing charge frustration. 
This frustration will spread over the whole lattice and 
it leads to the collapse of the long-range charge order. 
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FIG. 2: (Color online) (a) Spin structure factor ^(q) for vari- 
ous spin orders as a function of y, at Jaf=0.1 and A=1.4. (b) 
Typical MC snapshot of the spin configuration showing the 
short range FM domains (j/=0.094, Jaf=0.1, A=1.4), where 
the substitution ions are represented by filled circles. This is 
in the "FM cluster" regime (R2) discussed in the text, (c-e) 
Typical MC snapshot of the charge distribution. Here the 
circle area is proportional to the local charge density. High- 
and low-density sites (comparing with 0.5) are colored by blue 
and red, respectively, (c) The staggered CO pattern in the 
clean limit (Jaf=0.1, A=1.3). (d) CO pattern at j/=0.094, 
Jaf=0.1, and A=1.3. (e) CO pattern at j/=0.094, Jaf=0.1, 
and A=1.6, in the spin-glass regime (Rl). The locations of 
the substituting ions are the same in (b), (d) and (e). (f) 
Energy difference between the FM and CE ordered states as 
a function of A, for various values of y. 

Two typical MC snapshot of the charge redistribution 
are shown in Fig. 2(d) and (e) with intermediate and 
large A values. For the intermediate A case (here A=1.3), 
the charge density distribution is homogeneous except 
for some regions around the defects, while for the large 
A case (here A=1.6), the charge disproportionation is ob- 
vious although it occurs without a long-range ordered 
pattern. 

The intuitive idea described in the previous paragraph 
can be examined by means of a crude calculation: the 
phonons (classical distortions) are allowed to evolve freely 
in the MC sequence, while the i2g spin background is 
frozen into either the FM or CE patterns. The energy 
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FIG. 3; (Color online) (a)-(c) Calculated Cg electron density- 
of-states, DOS, for the two values of y indicated, at A=0.5, 
1.2, and 1.6. Jaf is fixed to 0.1. 

difference between the FM and CE phases, denoted by 
A defined as A=£'fm — Ece, is shown in Fig. 2(f). In 
the clean limit y—0, the energy difference A is relatively 
large, but it is rapidly suppressed upon increasing y. An 
appropriate choice of y and A may allow the energy of the 
FM phase to be even lower than that of the CE phase 
(A<0), and in these cases the ground state favors the 
FM order. This crude calculation illustrates the asym- 
metric impact of lattice defects on the stability of the two 
phases, suggesting a possible transition from the CE/CO 
phase to the FM phase. In fact, previous theoretical stud- 
ies on the A-site disorder effects have also confirmed the 
fragility of the half-doped CO phases .^>2^^>21i28^^^ 

Let us study now the electronic structure and trans- 
port properties of the system. The Cg electronic DOS 
provides insight on the effect of lattice defects. Figs. 3(a)- 
(c) show the calculated DOS at several values of y and 
A, for a fixed Jaf=0.1. The DOS at smaU A (A=0.5) is 
not qualitatively modified by the lattice defects, show- 
ing the anticipated robustness of the metallic state. For 
a large electron-phonon coupling A=1.6, the DOS shows 
a large energy gap at the Fermi level in the clean limit 
y=0, corresponding to the long range CE/CO phase. A 
substitution of ?/=0.094 clearly shrinks this gap, but still 
there are no states at the Fermi level (although there are 
some states close to it). However, for an intermediate 
coupling A=1.2, the most exotic features in the DOS are 
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FIG. 4: (Color online) (a) Jaf-A phase diagram at y=0.094 
compared with the clean-limit phase diagram Fig. 1 (a) (the 
blue lines are the phase boundaries in Fig. 1(a)). (b) Spin 
structure factor ^(q) at q=(0, 0) as a function of A, for the two 
substitution levels indicated, (c) Spin structure factor ^(q) 
at q = (0,0) as a function of Jap, for the two substitution 
levels indicated. 

obtained. In the clean limit, the gap is wide and obvious, 
but this gap completely vanishes at y=0.094, suggesting 
the stabilization of a finite DOS at the Fermi level and, 
if Anderson localization is not considered, metallic be- 
havior in the electronic transport. This lattice-defects- 
induced insulator to metal transition is similar to that 
induced by the A-site disorder M^^^^^^J^^^^S^ 

As a compact summary of the conclusions of this sec- 
tion, the MC calculated phase diagram at ?/=0.094 is 
shown in Fig. 4(a). Comparing with the phase diagram 
in the clean limit (see Fig. 1(a), and the blue lines in 
Fig. 4(a)), the FM metallic phase remains fairly stable 
and even expands slightly. This increase in the range 
of stability of the FM metal is important to rationalize 
the experimental results of Ref. 20. An interesting fea- 
ture of the phase diagram Fig. 4(a) is that the CE/CO 
regime fully vanishes, with the original boundary with 
the FM regime shrinking backward slightly. The orig- 
inal clean-limit CE/CO phase collapses into three sub- 
regimes: the FM metal and the two regions denoted in 
the figure by Rl and R2. Here, Rl (large Jaf and A) 
corresponds to a "spin glass regime" with short-range 
charge order and no visible FM order parameter S{Q, 0). 



This regime was described in previous publications, such 
as Ref. [23. In the more novel regime R2, the microstruc- 
ture consists of short-range FM clusters with various ori- 
entations (Fig. 2(b)). The charge order is suppressed in 
this regime (Fig. 2(c)), and this region is here called the 
"FM cluster regime". With increasing Jaf and A, these 
FM clusters/domains will be separated into even smaller 
domains, and eventually into a spin glass state. For this 
regime, there is relatively strong FM tendency, and a 
visible drop in 5(0, 0) (by varying A or Jaf as shown in 
Figs. 4(b) and (c)) is observed when the parameters (A 
and Jaf) cross the phase boundary between FM metal- 
lic and FM cluster regimes. In addition to their influence 
on the modification of the phase diagram, the lattice de- 
fects also smear the phase boundaries, implying inhomo- 
geneous tendencies in the FM cluster state. As shown 
in Figs. 4(b) and (c), the FM order parameter S'(0,0) 
decreases slowly with increasing A or Jaf, in contrast to 
the abrupt drop characteristic of the first-order FM-CE 
transition in the clean limit. Summarizing, the B-site 
lattice defect disorder replaces the clean-limit CE/CO 
phase by three different regimes: (i) A simple extension 
in parameter range of the competing FM metallic phase. 
In Figs. 4(b,c), this regime is between the A or Jaf where 
the original first-order jump from FM to CE in the clean 
limit occurs, to the clearly visible change in the slope 
of the 5(0, 0) curve (cusp) with further increasing cou- 
plings, (ii) The next regime is the FM cluster state (R2), 
already described, (iii) The following is the spin-glass 
state (Rl), also described before in detail in this section. 
Thus, we predict that half-doped manganites in the CE 
state could be destabilized in three different manners by 
B-site defects depending on how close they are to the FM 
metallic state in the clean-limit phase diagram. Qualita- 
tively, these results are similar to the effects caused by 
A-site disorder, and this similarity is a consequence of 
the previously unveiled fragility of the CE phasei^ 



C. Effect of electronic-density variations 

In the previous section, we have only considered the 
effects of the lattice defects introduced by the non- 
magnetic substitutions, which gave rise to results sim- 
ilar to those of previous studies that focused on A- 
site disorder. Now let us incorporate the effect of the 
eg electron-density variation due to the B-site substitu- 
tion. For the pure system, n>0.5 usually corresponds 
to FM order, while n<0.5 corresponds to the C-type 
AFM order, as shown in Fig. 1(b). After the substitu- 
tion, this density- variation-driven phase transition is still 
relevant i^i^ If the substituting ions have charge higher 
than -1-3.5 (ns>3.5), corresponding to 7ii->0.5, the spins 
present a FM tendencyi^i^ However, if ns<3.5, the situ- 
ation becomes much more complex since rtr<0.5 usually 
leads to the C-type AFM order instead of the FM one. 
For example, if the impurity is Cr/Al/Ga with charge 
-1-3 (ns=3), rir becomes less than 0.5 and the effective 
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electronic density will decrease to (0.5 — y) /(I — y) with 
increasing y. For 2/=3.1%, 6.3%, 9.4%, and 12.5%, the Cg 
electron-density n-^ in i?o.5^o.5Mni_yi3y03 drops down 
to 0.484, 0.467, 0.448 and 0.429, respectively. Therefore, 
the reduction of electronic density will compete with the 
FM tendency induced by the lattice defects discussed in 
the previous section. 

Let us recalculate 5(q) for different orders as a func- 
tion of y, at Jaf==0.1 and A=1.35, as shown in Fig. 5(a). 
Remarkably, a significant FM tendency is still present 
and the maximum 5(0, 0) appears at y~0.094, which can 
be considered as the optimal substitution level for the FM 
order. In contrast to the much reduced C-type AFM or- 
der in Fig. 2(a), here S'(0,7r) is partially sustained due 
to the charge density reduction, while on the other hand 
S{-k/2, 7r/2) vanishes quickly, suggesting the rapid disap- 
pearance of the CE order. A typical MC snapshot of 
the spin configuration considering the Cg charge-density 
reduction effect is shown in Fig. 5(c), where both the 
FM tendency and C-type AFM tendency are observed 
simultaneously, namely there are pairs of FM spin chains 
coupled antiferromagnetically as in the C state, and also 
small pockets of ferromagnetism. More specifically for 
the FM component, the calculated FM fraction in the 
spin structure factor in our simulations is about 17.5%, 
which is consistent with the experimental low-T magne- 
tization of Pro. 5Cao.5Mno.975Alo.o2503: 0.7/zb//.u., cor- 
responding to ~20% of the saturated magnetization.^*^ It 
should be mentioned that the above predicted optimal 
value of y is higher than the experimentally identified 
one. This disagreement may be ascribed to the A-site 
disorder already existing in Pro.sCao.sMnOa and other 
defects which are not considered in the present model. 
Although our results do not seem quantitatively accu- 
rate in this respect, we are confident that the qualitative 
tendencies have been captured in our calculations. 

As a conclusion of this section, the MC calculated 
phase diagram at y=0.094 with +3 impurity cations has 
also been calculated, and it is shown in Fig. 5 (d). Com- 
paring with Fig. 4(a), the effect of the electron density 
reduction is clear: the FM metallic regime shrinks while 
the FM cluster regime is relatively enlarged. For the 
regime above the dashed dot line, including portions of 
Rl and most of the R2 regime, the C-type AFM com- 
ponent remains robust (see 5(0, tt) in Fig. 5(b)). In 
addition, the FM order 5(0,0) in Fig. 5(b) decreases 
more smoothly with increasing A compared with the two 
curves in Fig. 4(b), indicating that the lattice defects 
and charge density reduction will both smear the phase 
boundary. Thus, the reduction of Cg electron-density sup- 
presses both the long-range CE and FM spin orders and 
enhances the importance of the C-AFM order. Therefore, 
the total combined effect of lattice defects and electron 
density reduction over the clean-limit CE state leads to 
an inhomogeneous state with coexistence of short-range 
FM and C-AFM spin ordering. 
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FIG. 5: (Color online) (a) Spin structure factor S'(q) for var- 
ious spin orders as a function of y (A=1.35) after taking into 
account the eg charge-density reduction effect, (b) Spin struc- 
ture factor ^(q) for two spin orders as a function of A. (c) 
Typical MC snapshot of the spin configuration considering 
the eg charge-density reduction effect at Jaf=0.1 and A=1.35, 
where the substitution ions are represented by filled circles, 
(d) Jap-A phase diagram at y=0.094 with nr=0.448. The 
dashed dot-line divides the phase diagram into two regimes, 
with the region of stronger A and Jaf having a stronger C- 
type AFM signal. 

IV. CONCLUSIONS 



Our Monte Carlo investigations reported here have 
shown that the role of the non-magnetic B-site substi- 
tution in manganites is rather complex. There are two 
main tendencies that compete: (1) The B-site substitu- 
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tion introduces lattice defects that break the CE-zigzag 
chains and causes charge frustration; (2) It also varies 
the Cg electron-density that leads to a density- variation- 
driven phase transition. In principle, both of these two 
roles are absent in the A-site disorder case. However, 
the lattice defects appear to induce similar results as the 
case of A-site disorder: part of the original CE phase 
regime in the clean-limit phase diagram is taken over by 
the FM metallic phase, or short-range EM clusters. This 
can be understood in the context of the previously dis- 
cussed "fragility" of the CE phase,— as compared with 
the robustness of the EM order. However, the concomi- 
tant reduction of Cg electron density suppresses both the 
FM and CE spin orders leaving behind a large inhomoge- 
neous area that consists of coexisting short-range EM and 
C-AEM clusters. Furthermore, the competition between 
lattice defects and electron density reduction gives rise 
to an optimized substitution level for the FM tendency, 
qualitatively similar as found experimentally. 

Summarizing, here we have investigated extensively 
the CE/CO destabilization experimentally observed in 
half-doped manganites due to a small amount of B- 
site nonmagnetic substitution, by using the two-orbital 
double-exchange model. Our calculations have shown 



that the CE/CO phase can be easily destabilized by lat- 
tice defects, leading to a variety of interesting possible 
states that include the competing FM metal, or FM clus- 
ters with or without C-AFM regions (depending on the 
valence of the B-site substitution), or a spin glass state. 
In particular, the surprising FM tendency observed here 
driven by the nonmagnetic substitution into the CE state 
is consistent with several recent experimental results. 
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